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Abstract 

We theoretically study dynamic properties of thin ferromagnetic films in contact with normal 
metals. Moving magnetizations cause a flow of spins into adjacent conductors, which relax by spin 
flip, scatter back into the ferromagnet, or are absorbed by another ferromagnet. Relaxation of spins 
outside the moving magnetization enhances the overall damping of the magnetization dynamics 
in accordance with the Gilbert phenomenology. Transfer of spins between different ferromagnets 
by these nonequilibrium spin currents leads to a long-ranged dynamic exchange interaction and 
novel collective excitation modes. Our predictions agree well with recent ferromagnetic-resonance 
experiments on ultrathin magnetic films. 

PACS numbers: 75.40.Gb,76.50.+g,75.70.Cn,75.30.Et 
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I. INTRODUCTION 



The spin transfer in ferromagnet-normal-metal hybrids causes a number of exciting phe- 
nomena. Equihbrium spin currents across a thin normal-metal (N) layer separating two 
ferromagnetic (F) films explain the oscillating RKKY-type exchange coupling jjj as a func- 
tion of the spacer thickness. This nonlocal exchange interaction may stabilize an antiparallel 
equmbnum co„figura«o„ ,„ag„e* multilayers 2|. wh.ch display a g,a,. .nagnetoresis- 
tance pj. For thicker N spacers or tunneling barriers, this static ferromagnetic exchange 
vanishes, but the F layers can still be made to communicate by driving a dc current through 
the system. For example, depending on the direction of the current flow perpendicular to a 
layered F\N\F spin- valve structure, a torque can be exerted on the magnetizations. When 



one magnetization is fixed, this torque favors either para! 



lei or antiparallel configurations. 



Transport in the dynamic 



and may lead to a switching of the other magnetization 

regime of moving magnetization directions has attracted relatively little attention, however, 

and is the subject of this paper. 

We have recently launched the idea of spin pumping in weakly excited magnetic nanos- 

tructures. A moving ferromagnetic magnetization emits spins into adjacent conductors, 

exerting a relaxation torque and transferring an angular momentum out of the ferromagnet. 

The ejected spin angular momentum can either scatter back, relax in a nonmagnetic spacer, 

or be absorbed by a second ferromagnet. In the first case, the magnetization dynamics 

are not affected; in the second case, the magnetization motion is (nonlocally) damped by 

spin-flip scattering processes in the normal metal; and in the last case, the two ferromagnets 

become coupled by an exchange of itinerant spins, which may result in collective excitation 

modes. The spin-pumping concept has proven to be fruitful in the understanding of both 

n n 

the nonlocal damping mechanism |6j] and dynamic exchange interaction [_7] in hybrid F\N 
systems. The functionality of magnetic devices, such as magnetic random-access memories, 
is strongly affected by relaxation characteristics of the magnetic media. The nonlocal re- 
laxation and coupling can be large and even dominant compared with other mechanisms in 
ultrathin films, and can be engineered by composition, geometry, and magnetic configuration 
of hybrid F\N systems. We hope that our theory will be useful as a designer tool, helping 
to realize future nanoscale magnetic devices with improved speed, operation threshold bias, 
and power consumption. 
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II. ENHANCED DAMPING IN SINGLE F FILMS 



Consider an F\ N bilayer as in Fig. ^ Without a voltage bias, no spin or charge currents 
flow when the magnetization of the ferromagnet is constant in time. When the magnetization 



direction starts precessing (as, e.g., under the influence of an api 



ied magnetic field) , a spin 



current IP"™p is pumped out of the ferromagnet into the N layer p . When the ferromagnetic 
film is thicker than its transverse spin-coherence length Age, d > Age = n/lkp — kp\, kp 
being the spin-dependent Fermi wave vectors, this current depends on the complex-valued 



parameter g^^ = g]} + igj\ the interfacial "mixing" conductance j8|,|9[, by 

ir-' = ^(d'n..^-gV'^] . (1) 
An dt dt J ^ ^ 

Here the time-dependent order parameter of the ferromagnet is a unit vector m(t), assum- 
ing a monodomain magnet with a spatially uniform magnetization at all times. A detailed 
derivation of Eq. based on the scattering- matrix theory is given in Refs. Alterna- 
tively, this result can be derived in the framework of magnetoelectronic circuit theory 
using only energy and angular- momentum conservation The dc conductance matrix is 
defined by 

^'^'^'=E['^-n-Cn(C'j*] • (2) 
mn 

Here r,^„ i^mn) is the reflection coefficient off the F\N interface for spin-up (spin-down) 
electrons on the normal-metal side, which can be obtained by ab initio band-structure cal- 
cu,at.n. Q. defined . E, be .e„o— to. M,.„ .a.p„e. i.erface. 

as explained in Ref. jl^, i.e., in the case of (globally) diffuse transport, the spurious Sharvin 
contributions have to be stripped off the real-valued conventional (Landauer-Biittiker) re- 
sistances for a definite spin a: 

while the inverse mixing conductance has to be corrected for the normal-metal Sharvin 
contribution only: 

±^±--L (4) 

The tildes in Eqs. (0) and (@)) denote the renormalized conductances, which reduce to the 
bare values (j21) when the numbers of transverse channels in the normal metal, N^, and the 
ferromagnet, Np„, sufficiently exceed the contact conductances g'^"' . 



The total spin current, I^, across the F\N interface also has a backfiow contribution, I^'^'^'^, 
in addition to the pumped current, IP"™p^ see Fig. ^ The total spin transfer in the steady 
state, 

I, = If "'P - I^"^'^ , (5) 

is determined self-consistently by the spin accumulation close to the interface. A finite value 
of the spin current, I^, would indicate the presence of a spin-sink mechanism in the normal 
metal by, e.g., spin-orbit coupling, when the spin angular momentum is transferred from 
the electron system to the lattice. In the case of small spin relaxation, the dynamically 
created spin accumulation in the normal metal may serve as a spin-battery device jU]. In 
the opposite spin-flip scattering regime, the angular-momentum loss of the ferromagnet by 
Is results in a Gilbert-type damping of the magnetic precession Q]. 

The spin current out of the ferromagnet carries angular momentum perpendicular to the 
magnetization direction, corresponding to a torque r = — I^ on the ferromagnetic condensate 
0]. Disregarding interfacial spin flips, this torque is entirely transferred to the magnetization, 
which is described by a generalized Landau-Lifshitz-Gilbert (LLG) equation P-Iisl 

dm dm 7 , . 

— = -7m X He. + aom X — + ^I. , (6) 

where 7 is the absolute value of the gyromagnetic ratio, ao is the dimensionless intrinsic 
Gilbert-damping constant, Mg is the saturation magnetization of the ferromagnet, and V 
is its volume. Referring to Eq. we see that the real part of the mixing conductance 
contributes to the damping just like the intrinsic bulk constant ao which is thus smaller 
than the total Gilbert damping a = Oq + a', whereas the imaginary part of the mixing 
conductance contributes like an effective field. The additional damping a' is observable 
in, for example, ferroma gne tic-resonance (FMR) spectra. Although not necessarily so for 



ferromagnetic insulators 
a good approximation real 



It L the mixing conductance for intermetallic F\ N interfaces is to 
12j | and therefore 



a = K- — — -— . (7) 



K = \ corresponds to the perfect spin-sink model, when all pumped spins relax in the N 
layer and the backfiow I^'^^'^ vanishes, /t < 1 corresponds to a finite backfiow [lfll |. 

The spin-sink capacity of clean normal-metal layers in contact with a ferromagnetic film, 
as in Fig. [TJ is usually governed by spin-orbit scattering processes at impurities or defects. 



As the spin-flip probability e = Tei/tso (deflned in terms of the elastic scattering time Xgi 
and spin-orbit relaxation time Tso) rapidly increases with the atomic number Z, e oc 



we expect a larger spin-sink effect [and therefore k in Eq. ((7j)] for heavier metals 
(impurities as well as hosts). The extent of the hybridization of the conduction bands with 
p oi d orbitals also plays an important role. In particular, clean noble metals, Cu, Ag, and 
Au, with predominantly s conduction electrons are poor spin sinks with correspondingly 
small e, but Pd and Pt, whose conduction electrons have signiflcant d character, have a 
high spin-orbit scattering rates and are efficient spin sinks. We note that heavy or magnetic 
impurities can turn an otherwise poor spin sink into a good one. The hierarchy of the Gilbert- 
damping enhancement has been measured for normal-metal buffers of Cu, Ta, Pd, and Pt 
(in order of increasing damping) , in agreement with the aforementioned arguments 3| ■ 
In particular, the effect of spin pumping on the magnetization dynamics was shown to be 
negligible in the case of Cu, while Pt was a nearly perfect spin sink, resulting in a large 
increase of the Gilbert constant; this was also experimentally confirmed by Ingvarsson et al. 

To further investigate the spin-sink effect of the normal-metal buffer, we next consider 
a more complex system, consisting of a double N layer attached to the ferromagnet, see 
Fig. 121 An interesting situation arises when the layer Nl is a bad _spin sink, such as Cu, and 
N2 is a perfect spin sink, such as Pt Mizukami et al. 0, 22] experimentally studied 
the FMR line width in permalloy (Py)|Cu|Pt composites as a function of Cu {Nl) width 
L. Next to L, there are three relevant length scales in the problem: the Fermi wave length, 
\fi the elastic scattering mean free path, Aei, and the spin-diffusion length, Asd- In the 
coherent regime, L < Aj?, we do not expect large effects of a dusting Nl layer on the mixing 
conductance, which means that the damping enhancement ((7j) is approximately governed by 
9^F\N2 (with a quantum modulation when L ~ Xp)- If L > A^, the spin transport across the 
Nl spacer may be described by the diffusion equation [l^, provided that either the spacer 
or the interface is disordered. In this regime, the spin backflow can be partitioned between 
the ferromagnet and layer N2., see Fig. El The relevant effective mixing conductance qIq 
then has to account for scattering at both interfaces, F\N1 and N1\N2, as well as in the Nl 
spacer: 

nr ~ ^^Ti ^ -^A^i + ~ ' (S) 

9eS 9Fm 9N1\N2 
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where Rni is the resistance (per spin, in units of h/e^) of the Nl layer, ^Ari|Ar2 is the one-spin 
conductance of the N1\N2 interface, and we assumed that L < Asd- If in addition L < Aei, 
the bulk scattering, Rm, can be disregarded and the total resistance is given simply by 
the sum of the interfacial contributions. Xp thus sets the length scale for a sharp drop in 
the Gilbert-damping enhancement, followed by an algebraic decay with the effective mixing 



conductance (jH)) for L > Ai?. Such a damping drop was reported in Ref. 
hybrids. 

The regime L > Xp was studied in Refs. 



23 for Py|Cu(L)|Pt 



21 



where a smooth algebraic decay of the 



damping enhancement was measured for L < Asd followed by an exponential suppression 
for thicker Cu spacers, in excellent agreement with our theory 10]. Mizukami et al. 



22 



also offered an explanation for their measurements, using a phenomenological theory due to 



Silsbee et al. l2c 



The theoretically calculated damping profiles, a(L), as a function of L in 



Refs. and j2^ are barely distinguishable. In our opinion, the phenomenological approach 
employed in the latter has severe limitations, however. In particular, the separation between 



localized and conduction electron spins is not justified for itinerant ferromagnets 

like the transition metals. Furthermore, this separation leads to wrong results for insulating 
ferromagnets which, we believe, can generate a spin current into adjacent nonmagnetic 
conductors in the same way as a conducting ferromagnet llOll. in contrast to the prediction 
of Ref. Q- 

The good fit of the experimental a{L) in Ref. j22] reflects the diffuse nature of 
transport in the copper layer rather than the spin- injection mechanism at the F\N interface. 
The magnitude of the spin current as calculated in Ref. j22] is governed by a fitting parameter 
(r in their notation) which cannot be verified or predicted by ab initio calculations, unlike 
our spin-pumping expression (|T)). 

The generality of our approach is demonstrated in the following by showing its validity 
for a qualitatively different system of two magnetic layers separated by a normal-metal film 
with negligible spin flip. We find that the same mechanism (namely the spin pumping) 
is responsible for the damping enhancement in coupled ferromagnets as well as the F\N 
structures discussed so far. 
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III. DYNAMIC COUPLING IN SPIN VALVES 



In F\N\F spin valves, the spin pumping causes qualitatively different effects in addition 
to the ones just described. Consider a system shown in Fig. |21 In the following we take 
L > Xp, so that quantum coherence (in particular the static exchange coupling) can be 
disregarded, but L < Ad, so that the spins of electrons are transferred between the two 
ferromagnets ballistically. 

The total spin current pumped into the normal metal, 

has contributions from both Fi\N interfaces, see Eq. (Q). Here the small imaginary part 
of the mixing conductance 'l2] is again disregarded. In a collinear configuration, the spin 
accumulation fig induced by the spin pumping is perpendicular to the magnetizations. This 



simplifies the expression for the bias-driven spin current to 
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jback^^Ii+ili^^^ (10) 

A steady state is established when the two spin currents, Eqs. and ()1U|) . cancel each 
other: I^^'^'^ = Ip^i™p. The spin accumulation is then given by fig = 47r(^P + ^2^)~^IP"™p. 
In general, 1^2 = — Isi, by the conservation of angular momentum, but the net spin current 
Isj = Ig"™'' — I^J"^*^ across a single interface does not vanish. For the left (j = 1) interface we 
find 

The normal-metal layer in our model scrambles the incoming spin current and divides it back 
over both ferromagnets. [For very transparent interfaces the scrambling is only partial, but 
the present treatment is still adequate if we properly renormalize the conductance parameters 
as in Eqs. ^ and (jH).] The second ferromagnet can thus cause damping of the precession 
in the first magnetic film, and vice versa. A ferromagnetic layer can therefore serve as an 
efficient spin sink, just like a normal metal with strong spin-flip scattering. 

We first discuss the implications of the spin transfer (|1H) when the second ferromagnet, 
F2, is stationary, m2 = 0. In the FMR measurements, this is the case when one ferromagnet 
is in resonance, whereas the resonance frequency of the second ferromagnet is sufficiently 
different because of different magnetic anisotropies. The spin torque r = — I^i on the first 



ferromagnet {Fl ) then has the form of the Gilbert damping when added as a source term to 
the LLG equation ©• In this case, the dynamic couphng of the ferromagnetic layers simply 
leads to an enhancement of the Gilbert-damping parameter with respect to its intrinsic 
value, exactly like in the case of single F films. The damping enhancement, following from 
Eq. of the Fl layer therefore reads 

A7rM,VCgl^ + g^' 

(with 7, Ms, and of Fl). Eq. (IT^ satisfactorily explains the increased viscous damping 
observed in Fe|Au|Fe spin valves j2J|. The physical nature of the damping in single F films, 
Eq. ((Tj), versus spin valves, Eq. ()12|) . is quite clear. The angular momentum is first driven 
out of the ferromagnet via spin pumping. In the former case (single ferromagnet), the 
subsequent spin-orbit processes in the normal metal relax injected spins and, as a result, 
slow down the coherent motion of the ferromagnet. In the latter case (spin valve), the spin 
angular momentum is transferred from ferromagnet Fl and into F2. This spin current is 
absorbed by asynchronously driving the magnetization dynamics of F2. Therefore, both 
the normal-metal and ferromagnetic spin sinks act as an external brake to slow down the 
precession of the resonantly excited magnet. The spin-sink efficiency of the normal metal 
is characterized by the spin-flip probability e and its thickness L which is lumped into 
the parameter k in Eq. ((Tj). It follows from Eq. (fT^ that the efficiency of the adjacent 
ferromagnet F2, on the other hand, only depends on the F2\ N mixing conductance (as long 
as the second ferromagnet is thicker than its transverse spin-coherence length Age). 

In the following we consider dynamics of the coupled F1\N\F2 system when both mag- 
netizations are allowed to precess, i.e., when the ferromagnet resonance conditions are close 
to each other. By augmenting the LLG equation (jHI) for mi by the spin current source term 
and with small variables u,; = m.j — hi, where \ui\ ^ 1 and Uj ± hj, we obtain the 
linearized expression (assuming a circular precession) 

ill = ^^ihi X ui + aohi x lii + a[ (hi x lii — h2 x U2) , (13) 

where cuj is the resonance frequency of ferromagnet Fi, hj is the unit vector in the direction 
of its effective magnetic field, and a- = •jihgl^gl^ /[iirMsiVi^gl^ + gl^)], subscript i labeling 
corresponding quantities of Fi. The dynamics of the magnetization direction 1112 is obtained 
by exchanging subscripts 1 ^ 2 in Eq. (|T^ . 
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For a spin valve in the parallel configuration, hj = h, and identical resonance frequen- 
cies, Ui = uj, the quantity u = ui/a[ + ^2/0.2 (which, up to a scaling factor, is the rf 
component of the total angular momentum) is affected only by the intrinsic bulk damp- 
ing, li = uj\i X u + agh X u, while the difference Au = ui — U2 relaxes according to 
Au = tuh X Au + ah x Aii with the enhanced (Gilbert) damping constant a = + a'l + ot'2- 
The dynamic coupling in the antiparallel configuration as well as in the parallel configura- 
tion when the resonance frequencies have a large mismatch, Acj ^ (ao + ct'i + «2)'^) favors 
relaxation of each magnetization towards the equilibrium configuration. In this case, the 
two modes corresponding to excitations of either ferromagnet decouple, each having an en- 
hanced damping parameter a-i = + a[. This picture explains in detail the FMR profiles 
measured on the Fe|Au|Fe spin valves, not only in both symmetric and very asymmetric 
limits discussed herein but also in the intermediate regime of closely matched but different 
resonance frequencies in the parallel alignment 0. 



IV. CONCLUSIONS 



We presented a simple yet (for diffuse systems) complete theory of spin-current emission 
into adjacent conductors, which is generated by the time-dependent ferromagnetic order 
parameters. In the case of a single magnetic film in contact with a normal-metal layer, spin- 
flip processes in the latter viscously slow down the coherent precession of the ferromagnet, in 
accordance with the Gilbert phenomenology. When a magnetic layer is connected to other 
ferromagnets via Ohmic contacts, the latter can act as a nonlocal brake by relaxing the 
transverse angular momentum of the incident spins. Our framework formulated in terms 
of spin pumps and spin sinks successfully explains ferromagnetic resonance experiments on 
various multilayer systems. 

In addition, this approach predicts a dynamic exchange coupling between ferromagnetic 
films connected by normal metals. Precessing magnetizations feel each other through a 
metallic spacer by exchanging nonequilibrium spin currents. Components of the ferromag- 
netic magnetizations that precess incoherently (i.e., out of phase) experience an additional 
damping obeying the Gilbert phenomenology. In symmetric spin valves and periodic super- 
lattices, the dynamic exchange coupling causes a weakly damped, synchronized collective 
motion of the magnetic order parameters. 
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FIG. 1: Schematic view of the F\N bilayer. Precession of the magnetization direction m(t) of the 
ferromagnet F pumps spins into the adjacent normal- metal layer N by inducing a spin current 
jpump_ xhis leads to a buildup of the normal-metal spin accumulation which either relaxes by 
spin- flip scattering or flows back into the ferromagnet as I^^^*^. The N layer here is not an ideal 
reservoir but rather a film of the same cross section as the magnetic layer F; the spin accumulation 
is position (x)-dependent. 



FIG. 2: Same as Fig. ^ but now the normal-metal system is composed of a bilayer N1\N2. Fer- 
romagnetic precession pumps spins into the first normal- metal layer Nl. The spin buildup in Nl 
may flow back into the ferromagnet F as spin current I^f'^'^, relax in Nl, or flow into the second 
normal-metal layer N2 as spin current I^f'^'^. The spin accumulation in N2 is disregarded since the 
layer is assumed to be a perfect spin sink. 



FIG. 3: Two ferromagnetic films Fl and F2 connected by a normal-metal spacer N of width L. 
When the magnetization directions mi and m2 precess, spin currents I^"™'' and 1^2™^ pumped 
into the normal metal. The spin-dependent chemical-potential imbalance in N causes the backflow 
of spin currents t^f^^ and I^f'^'^. 
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Figure 1. Y. Tserkovnyak et a/., J. Appl. Phys. 
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Figure 2. Y. Tserkovnyak et a/., J. Appl. Phys. 
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Figure 3. Y. Tserkovnyak et a/., J. Appl. Phys. 
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